home *** CD-ROM | disk | FTP | other *** search
/ WINMX Assorted Textfiles / Ebooks.tar / Text - Mathematics - Numerical Mathematics and Computing (F).zip / mctst1.f < prev    next >
Text File  |  2002-06-11  |  708b  |  28 lines
  1. C
  2. C PAGE 347: NUMERICAL MATHEMATICS AND COMPUTING, CHENEY/KINCAID, 1985
  3. C
  4. C FILE: MCTST1.FOR
  5. C
  6. C VOLUME OF A COMPLICATED REGION IN THREE-SPACE BY MONTE CARLO (RANDOM)
  7. C
  8.       DATA  L,M,N/256,0,1000/ 
  9.       DO 2 I = 1,N
  10.         X = RANDOM(L)       
  11.         Y = RANDOM(L)       
  12.         Z = RANDOM(L)       
  13.         IF((X*X+SIN(Y) .LE. Z) .AND. (X-Z+EXP(Y) .LE. 1.0))  M = M + 1
  14.         IF(MOD(I,100) .EQ. 0)  THEN   
  15.           VOL = REAL(M)/REAL(I)       
  16.           PRINT 3,I,VOL     
  17.         END IF    
  18.    2  CONTINUE    
  19.    3  FORMAT(5X,I5,F10.5)   
  20.       STOP
  21.       END 
  22.   
  23.       FUNCTION RANDOM(L)
  24.       L = MOD(16807*L,2147483647)
  25.       RANDOM = REAL(L)*4.6566128752458E-10      
  26.       RETURN
  27.       END 
  28.